
STOP 



Early Journal Content on JSTOR, Free to Anyone in the World 

This article is one of nearly 500,000 scholarly works digitized and made freely available to everyone in 
the world by JSTOR. 

Known as the Early Journal Content, this set of works include research articles, news, letters, and other 
writings published in more than 200 of the oldest leading academic journals. The works date from the 
mid-seventeenth to the early twentieth centuries. 

We encourage people to read and share the Early Journal Content openly and to tell others that this 
resource exists. People may post this content online or redistribute in any way for non-commercial 
purposes. 

Read more about Early Journal Content at http://about.jstor.org/participate-jstor/individuals/early- 
journal-content . 



JSTOR is a digital library of academic journals, books, and primary source objects. JSTOR helps people 
discover, use, and build upon a wide range of content through a powerful research and teaching 
platform, and preserves this content for future generations. JSTOR is part of ITHAKA, a not-for-profit 
organization that also includes Ithaka S+R and Portico. For more information about JSTOR, please 
contact support@jstor.org. 
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AND ANDREWS' PROFESSOR OF ASTRONOMY IN THE UNIVERSITY OF DUBLIN. 



Read April 21, 1817. 

THE difficulties that occur relative to the investigation of the mean 
motion of the apsids of the lunar orbit are well known. Two cir- 
cumstances have principally led me to offer to the Academy the 
following paper on this subject. It affords me an opportunity of 
introducing a peculiar method of integration which I hope soon to 
illustrate more fully elsewhere. This method is I think deserving 
of the attention of mathematicians, and as applied to the integra- 
tion of the equation of the lunar orbit affords very convenient results. 
The other circumstance alluded to, was a desire to improve a parti- 
cular step that appeared defective in the lunar theory as given by 
several authors, and recently by M. Laplace. These circumstances 
and some deductions that offer themselves in the following investi- 
gation will be best understood by a brief detail relative to this subject. 
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The computation of the mean motion of the lunar apsids is by 
far the most important point in which the Newtonian Theory of the 
Moon, as given in the Principia, appears defective. Newton him- 
self, in the latter editions of the Principia, seems to have aban- 
doned the attempt to reconcile or rather deduce from Theory the 
motion given by Observation. 

In the first edition he had made the attempt, after stating his re- 
sults (Scholium, p. 462.) ; he adds, " computationes autem, ut nimis 
" perplexas & approximationibus impeditas, neque satis accuratas, 
" apponere non lubet." 

It may be presumed that he found his method on examination 
inaccurate, otherwise it cannot be doubted he would have no- 
ticed it in the subsequent editions, and given, if not the method 
in detail, the results. 

Machin appears the first after Newton who attempted this 
problem ; the inadequacy of his solution, and of those of some 
ethers of the same nature, will be noticed further on. 

Clairaut, in 1748, had the honour of giving the first exact 
solution according to the principles of the Newtonian Theory of 
Gravity, after he had, in the Memoirs of the Royal Academy 
of Sciences at Paris, announced that the Newtonian law was in- 
exact, inasmuch that the mean motion of the lunar apsids deduced 
from that law did not agree with observation. 

Clairaut's result was confirmed by Euler, D'Alembert and 
Mayer, and subsequently by other mathematicians. Their re- 
searches, however, being directed more towards a general theory 
of the lunar motions, than towards the particular question of the 
mean motion of the apsids, are so complicated, that the exact thread 
of reasoning respecting this motion cannot without difficulty be 
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traced. The integration of the principal equation is usually com- 
menced by supposing a formula which depends on the knowledge 
of that integration or on the result of observations. Thus in 
the Theory of the Moon by M. Laplace, he has, * u being the 
reciprocal of the Moon's Distance from the Earth, 

u = A'(Uy') { J + x * y* + e cos ('—«•) &c -^ 
and remarks 

" Cette valeur de u suppose l'ellipse lunaire immobile ; mais 
on verra bientot qu'en vertu de Paction du soleil, les nceuds St le 
perigee de celte ellipse sont en mouvement. Alors, en designant 
par (1 — c) » le mouvement direct du perigee, &c. 

This is his first approximation, but certainly the first approxima- 
tion should be the former, and the second should be deduced 
therefrom by a regular process. The result undoubtedly con- 
firms this hypothesis, but it seems more consonant to the usual 
steps of mathematical reasoning to deduce one from the other, 
this is an object in the following investigation, in which also 
the mean motion of the perigee is computed by confining the 
process principally to this point, and therefore will be easily 
intelligible to those who may be unwilling to encounter the 
formidable calculations necessary for the complete lunar theory. 

The method of integration which is applied to the differen- 
tial equation of the lunar orbit is peculiarly convenient for 
the above purposes. With respect to this method it does not seem 
necessary here to remark more respecting it, than only to men- 
tion, that it principally derives its convenience from certain theo- 
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* Mec. eel. p. 187. Liv. 7. Tom. 3. 
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rems for finding Fluxions per Saltum. The method will be 
easily understood without the theorems. These and important 
applications thereof it is intended soon to give in a separate 
volume. 

In order to simplify the computation as much as possible, the 
plane of the lunar orbit has been supposed coincident with the 
plane of the ecliptic, the orbit of the earth without excentricity, 
and the approximation has been only carried to the first power 
of the excentricity of the lunar orbit. These circumstances have 
little effect on the quantity of the mean motion of the Lunar 
Perigee. 

The quantity of the motion found is expressed in terms of 
the quotient of the periodic time of the moon by the periodic time 
of the earth, and thereby are satisfactorily shewn the erroneous 
conclusions of Machin, Walmsley, Frisi and Matthew Stewart, 
who imagined that the mean motion of the apsids could be in- 
vestigated by considering the moon acted on only by a centripetal 
force, the mean tangential force being = o. This is of some im- 
portance, as authors have recently referred to these solutions as 
exact. Professor Playfair, indeed, in his outlines of Natural Philo- 
sophy, published in 1814, speaks (vol. 2. p. 26 1) with some doubt 
on the subject. After giving Dr. Stewart's result, and referring to 
those of the others, he says, " The result of these investigations, 
" therefore, agrees nearly with observation, but it cannot be 
** denied that the principle on which they are founded is liable 
" to some objections, so that if it were not for the information 
" derived from the direct solution of the problem of the three 
" bodies, it might still be doubted, whether the principle of gravity 
" accounted exactly for the motion of the Moon's apsids." 
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It is entirely by accident that their results are exact in the case 
of the Moon. Had the periodic time of the moon been different 
from what it is, observations would have pointed out the error of 
their conclusions. 

The differential equations and the expressions for the forces acting 
on the moon might have been taken from Laplace or other authors, 
but it is hoped that the manner in which the differential equations 
are deduced, may afford an excuse for the introduction of the in- 
vestigation. The method of limits, the undisputed invention of 
Newton, seems deserving of more attention than is nowpaid to it. 

The differential equations are at once deduced without the inter- 
vention of rectangular co-ordinates which in the present case would 
unnecessarily have added to the length of the process, but in other 
respects advantage has been taken of the facilities offered by La- 
place's investigation of the Lunar motions which may be consi- 
dered as the most satisfactory and as the most accurate that has 
appeared. 
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(I) 



Investigation of the differential equations of the orbit described 
by a body moving in a fixed plane about a fixed centre, when 
acted on by two forces, one directed to that centre, and the 
other in a direction perpendicular to the radius vector. 



Let APQ be the orbit described 
about the centre E; let PE =zx ; AEP 
= » ; velocity at P—v ; the force at P in 
direction PN perpendicular to PE=P ; the 
centripetal force at P=R ; area AEP=s 
and the time of describing PQ = At; the 
sign a expressing a difference. 

Then if PT be drawn perpendicular to q 
EQ we may thus proceed to investigate 
the relation of v and x 



_ = hnat ___ ~ hm. ^ = hm. ; 7 == FF 

P Q d A P 

now v = lim. — = -jj- or dAP = vdt 




= lim. 



\ PT* 



Therefore lta« - lim. ^ = ^x 
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(2) 

To investigate the value of ~ we may consider z a function of t, 
then by Taylor's Theorem 

(3) 



*• = **«+&*«•+&**+•* 
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As may also be considered as consisting of two parts, one 

part (= j t A t ) is independent of the two forces R and P, which 

between P and Q are variable both in quantity and in direction, 
and the other part arises from the action of these forces during the 
time a t. This part we may suppose also produced by an equivalent 
constant force acting during the time a t in a direction perpendi- 
cular to PE, and we may represent this force by P + q a t, q 
being a variable but determinate quantity, then by the property 
of uniformly, accelerated motion 

a* = ^a*+^-(P+<?a*)a* 2 (4) 

comparing equations (3) and (4) 

^+£U« + &<>•=§ ( P +9*<> 
or taking the limits of this equation 

d*z Px J9 Pxdt 1 .„ 

jTi = T or d*z = - r (5) 

But^ = lim.f->^ord S = ^ (6) 

Multiplying together equations (5) and (6) 
dzd 2 z = ~ x^dvdt* 

4 

and by integration 

dz* =(h+f?£ d,)dl* (7) 

h being a constant quantity. 

Then from equations (7) (2) and (i) we obtain 

7 d X 

d , = -, (8) 






4 h+f$Px*df 

also by equations (7) and (6) 

at zz : . — =■-==■ (a) 
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But it is necessary to eliminate v from equation 8. 
Now ~ = Km. ~ = Force in direction of the Tangent 

at A * 

ni . PT t> r QT _ D dt t> dx 

= P hm. -s^ — Jl hm. -sa = P a? t-j-b — It 



pQ "•""""jPQ "-"dAP " dAP 

AT 

dt 



Therefore because v = —r- 



vdv = Pxdt — Rdx 
or by integration 

t>* = / 2Pxd» — f2Rdx + k t k being a constant quantity Qj} 
Hence equation (8) becomes 

di> d t 

x \ l ULflEzl. Hub* u_ i_ f &") 

This equation of first fluxions is not so convenient for integrating 
by approximation as the equation of second fluxions which can 
easily be deduced from it, as follows. Both, however, will be 
used in the subsequent investigations. 

Substituting u = — and making y = fPu~ d* +fRu "~ d u + ~ k 
and it) = 2 h + f P w"~ 3 d », we obtain by squaring, &c. 

**(t- «»)=*«» (9) 

differencing this equation, making d v constant 

(10) 



dv y di» d 1 u 



2 terfu — 2w" au "~ di % 

but by equat. (9) 

yd™ _ JL x JLOL — (iaL + m 2 ^ — 

id 1 du to wdu \d>* J ndu 

therefore by equation (10). 

£(*-<£ + -)£) — = £ <"> 
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Now! " *»~i+ fiM ~* (11) 

anc j sou* \ aw __ Pu d u Pu d > 

\d?~l ru *)du d» "cHT 

hence we easily deduce from equation (11) 

The equation a is equivalent to the first of the equations (L) of 
Laplace,* and the equation (c) is equivalent to the second of the 
equations (L) when s = o. 

(II) 
Application to the Lunar Orbit. 

1. If we suppose the moon only attracted towards the earth, 
P = o, and R varies inversely as the square of the distance from the 
earth. Let R = ^., M representing the sum of the masses of the 

moon and earth. Let also a(l+e)anda(l — e) represent the 
greatest and least distances of the moon from the earth on this hy- 
pothesis : these distances must be invariable, because the centripetal 
forces (the force P being = o on this hypothesis) are equal at equal 
distances on each side of the apsids. We can obtain the values of the 
constant quantities h and k of the equation (6) in the following 

manner. Equation (9) becomes v l — — / — f- ^ = — + k 

therefore V and V representing the velocities at Perigee and Apo- 
gee respectively 



When P=:o equation (7) gives 
vol. xin. G 

* Mec. eel. p. 181. torn. 3. 
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*1) = h, therefore at the Apsids 

4dz l \h 



V' = 



dt * X «/«*•* ~" <*»* 2 
2l » 4/« 



hence 2m + k - 



a(l— e) T ~" rt « (!—(•)' 



2M + ft =- 4A 



These two equations give 
h = ; o (l— e 2 ) M 



M 



a 



Substituting these values in equation (6) and making P = o 
andfl = 4 






i!.i (15) 



Va; a 



I 



The integration of this equation (making the longitude of Peri- 

gee = *■), gives » = arc (cos = — — 7 ) + «• 

/ x o(l— e l ) 1 
or cos f» — it) = - —- - 

* ' ex e 

and therefore 

s _ g=. e Ii— : (16) 

This is the equation of an ellipse, the centre of force being in the 
focus. 

The second law of Kepler follows from this conclusion. 

It may be remarked here that the first law of Kepler follows 



from equation (7) when P = o, for it gives z — ht, as we may sup- 
pose z and t to commence together, and therefore the areas 
about a fixed centre are proportional to the times of describiug 
them. Also by help of equation (16) we deduce a conclusion that 
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will be of use hereafter, and which also proves the third law of 
Kepler. 

Because t = ;-= gl— -— ; consequently, 

since by equat. (16) the orbit is an ellipse, 

the periodic time is as M g gffi (17) 

and hence about the same centre of force, neglecting the masses 
of the revolving bodies, the squares of the periodic times are as 
the cubes of the greater axes. 

2. Supposing the moon also acted on by the sun. A small force 
P exists, and also a small alteration of the force R takes place, and 
then the integral of the equation (Z>) is less easily obtained 
than that of equation (c). Now as the values of P and JR de- 
pend partly on the relation of w and t, we cannot exactly ascer- 
tain the former without knowing the latter ; therefore previously to 
the integration of equat. (c), we must use approximate values of 
these quantities obtained by the integration of equation (b), when 
P = o and R — -f 
By equation (16) 

l+e cos (» — w) 
M = a(l_e*") 

or if we regard only the first power of the excentricity, 

« = a (1 +c cos (»—*)) ( 18 ) 

also h = \aM 

Therefore when P = o, equation (a) gives 

4 -. 

dt = 1± = a -4(\-2 e cos O-O) 

(l+ec<«(»— *r)) ! M* itf' 

g2 
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By integration, 



i 
t =°— (» — 2<? sin (v— t)) wliich requires no correction, as the 

origin of t is arbitrary. 

The above is relative to the moon about the earth ; but if we 

suppose v and a to be the same quantities relative to the earth 
about the sun, or rather the sun about the earth as » and a are 
relative to the moon about the earth. 

supposing the earth's orbit without excentricity and S (the mass 
of the sun) very great in comparison of M. 





3 

% 








3 


Hence ^ 


{v-%e 


sin {v — sr)) = 


a » 

X. 




iH« 








S* 


or 


making 


m — 


i I 
a* S 1 


Periodic time moon 
Periodic time earth 


j 


= m v — 


- 2 em sin (v- 


-c) 





nearly (19) 



Having thus obtained approximate values of u and ' in terms 
of v y we proceed to find the approximate values of P and R 
for substitution in the equation (c) 

In the preceding figure E, P and S represent the places of the 
earth, moon and sun. 

SE = d; EP ==-lau SEP = v—{ 

The attraction between S and E will be represented by j-^> tne 

mass of the earth being neglected in comparison of that of the 
sun. 
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The difference of the forces of S on P and on E in a direction 
parallel to SE = S (§,-'—) 
The effect of this difference in direction EP opposite to JEa: 

S (U — hO cos O-'O 

and in direction perpendicular to PE = S fjp- s — uTi) sin O—'O 

PE 

also the force of S in direction PJ3 = S x jrr; 

Hence R^^+S (^) — S (|^ — i_) cos (»-.;) 

and P = 8 (gj — Lg;) sin (,—',) 

or in consequence of the substitution m = - T — i? if we represent 
M the sum of the masses of the earth and moon by unity 

t, 1 , m*a* <PE ,SE 1 \ /- /x ? 

11 = ?F+ "ST- [spi ~ (sF ~ s&) cos 0—0 J 

Pm l a* /SE 1 \ . ^ / x 
= — (ST» - JT*) sm 0-0 

now sp =: V a ~^ - 7 ( cos ^ = « y/^-^^^b 

hence neglecting quantities of the order j 

&P 3 — ^ V 1 «« ./ 

Therefore if we neglect quantities of the order ^~ 

R = M * + £L__ ?f. a cos 2 (>-»') = m 2 — fr (1 +3 cos(?»— 20) 

P = — -T- cos (*— ») sm (>- J) = — 9— sm (2*— 2») 

„ fPer. Time Moon\ • 1 i 

now as m* = ( Pgr . ^ Earth ) 2 nearly = — nearly 

, , . 1 Parallax Sun i 1 i 1 _„_ i 

and a being unity j = Pu)alkx Mo - nearly = m nearly ; - may be 
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considered of the order m 2 . Therefore in the above values of It 

aud P quantities of the order m* are neglected. 

These values of R and P are next to be substituted in equation (c) 

*• 5? = - Th^? mn (2 " - 2") = 33- «n (2,-2,) x 

(i — 4 e cos (V— *)) by equat. (18) 

But sin (2, — 2 ,) = sin (2 » — 2m v + 4 em sin (v— sr)) 

= sin (2 x— 2 mv) + 4em sin (y — x) cos (2 v— 2 m v) 

— sin (2 v— 2m ,) +2 em sin (3 » -2 m»— *) — 2em sin (,— 2m, + sr) 




sin (2 v — 2 m ,) 

— (2 — 2m) e si 

— (2 + 2m) e sin (,— 2 m * + «•) 



hence ^ = — 3 -g-<( - (2-2m) e sin (3 ,-2m ,-*) 



consequently 



2 

cos (2, — 2 my) 



2— 2m 
I r? Pd, 3»; l « 1 2— 2»« ,„ a A 

2TJ "ST = -5T «{ ~ 5=^ e cos (3, -2m, -V) 

— f=S ccosO— 2m» + *) 

Without considering the disturbing force of the Sun, we found 
M being unity that h zz •* « (1 — e 2 ) = \ a neglecting the se- 
cond power of the excentricity. If therefore we make h = ~a + hf 
h' will be a quantity of the order of the disturbing force or of 
m 2 

Also without the disturbing force of the sun, equation (c) 
becomes 
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2^- cos (2 v- 2m 

, /d*u \ siPdt Sm J 1 /9 2>»\ 

hence (- + „)/_=_-<> _(§-|?) e cos (3 p-8m.-») 

—('££) e cos (" 2W+*) 

neglecting quantities multiplied by ti m 2 as being of the order m 4 

2. ££~ = 5 m* sin (2 ,-2 >>in (,-x) (1-3* cos (_r)) 

3m « Te cos (v— 2m >-+*•) 
— ilj — e cos (3 v — 2mv — t) 

3. Lastly - ^ = - ^, + - - — e cos (v— * ) 



f cos 2 >-2m» 



+ 3 _^1 J — —^ e cos (»-2!»»+*0 
£ __ 3^ e cog ^ 9 __ 2 m»-y) 

Substituting these quantities in equation (c) it becomes 

r-r^u r—J7 + s— tt- ecos(V— sr)+ 5—1 1 lcos(2»— 2m») 

14-2w-. j- I— m ~ 

9 " * 2 jL( e cos(,-2m, + *-)-^j t_Ae cm (3,-3m>) 

= o. (d) 

Each of the terms after the four first terms are contained in the 
form m*K cos (en — r) : it is therefore required to integrate an equa- 
tion of the form 

— J,, u — G + m? K (cos c»— t) + #c =o 
G and K being constant quantities. 
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(III.) 

Integration of the Equation. 

Jt + u — Gf -f- m?K cos (cv— v) = o 

G and K leing constant quantities y and 

1. c— any number greater or less than unity, 

2. c= 1 and m 2 of any magnitude. 

3. c = 1 and m 2 « very small quantity. 

The subsequent mode of integration will be more readily un- 
derstood by a short prefatory explanation.* 

Let us suppose we have in any manner arrived at an equation. 

d m <p (x, y)—o, when a; = o and dx is constant, m representing 
any number not less than w, we can then conclude 
<p {x, y)~ Ci+ c 2 x + c 3 x 2 — + c n x n ~ l 

The equation d m <p (x, y) = o may be called the mi' 4 particular 
fluxion of this equation with reference to x=o. 
If we have only equations of particular fluxions commencing with 
the order n, the values of c i} c 2 , c 5 , &c. are arbitrary. But if 
there be preceding equations of particular fluxions not contained 

* The mode of integration here shortly explained derives its advantage from the method 
of finding fluxions per saltum. Theorems for this purpose were given by me in a paper read 
before this Academy in the year 1798, and published in the Seventh Volume of the Trans- 
actions. These theorems were considerably extended, and applied both to the direct and in- 
verse reduction of analytical Junctions in a work which I had prepared for the press, but hav- 
ing adopted a notation differing both from the fluxional notation, and that of the differential 
calculus, I was for many years deterred from publishing it. Lately, however, I have again 
resumed the subject, changing the notation into the usual notation of the differential calculus, 
and I hope soon to offer the result to the notice of mathematicians. The method of integra- 
tion here used belongs to that division of the work entited " On the inverse reduction of 
analytical functions." 
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in the form d m <p (x, y) = o, such equations will serve lor determining 
certain of the values of c, , c 2 &c. 

It will be generally found more convenient in reference to the 
converse of theorems for finding fluxions per salt urn to use the 
fluxional equation. 

dm (p (x, y) 

1 . 2 . . m ' — ° 

which may be denoted by 

d m <p (x, y) = o, the denominator being understood by the 

line under d and by the order of the fluxion, so that 

d m <p {x, y)=o 
may with respect to the value of x = o and the denominator be cal- 
led a particular divided fluxional equation, of which the integral 
is 

<p (x, y) = Ci + c 2 x + c 3 x 2 c„ x H ~ l 

This mode of proceeding is applicable to the summation of many 
series, to finding generating functions, to the integration of equa- 
tions of finite differences and of many fluxional equations, more 
particularly those equations called linear equations. 

To proceed with the equation (e) 

It fluxion gives 

d 3 u + d u d r = m* Kc sin (c»- *) d ** (!) 

1. Taking the n—3 particular fluxion when *=o and d i is con- 
stant, n being even and not less than 4. 

n— 4 

d n u + d n - 2 udy' i = m 2 (—1) * Kc"--' cos *t di 1 (2) 

2. Taking the n— 3 particular fluxion when ^oanddr is con- 

VOL. XIII. H 



42 
stantj n being odd and not less than 3. 

n— 3 

d n u + d B ~ 2 tt d v 2 = — m 2 (—1 ) T" Kc"- 2 sin * d »" (3) 

Let a be a function of u 

such that when >=:o 
d n m = d" z 
Then equations (2) and (3) become 

n — 2 

d« g+dr~*z d i> 2 = (— 1 ) ~m 2 Kc n ~ 2 cos * d »" (4) 



d- 2 + d*-** d » 2 = -(— 1 )Tm* Kc"~* sin * d »■ (5) 

Now when v=o 

the nth particular divided fluxion of ^ , , 

when n is odd and not less than 3= d» 3 d n ~ 5 l zf^i 

— l)-i~ "•* and when n is even = o 
and the wth particular divided fluxion of 



v 3 when n is even and not less than 2 = ( — 1)2 c dv i 
i+c*? t-when n is odd = o J 



Hence equations (4) and (5) are deduced from the integral 
(l+»0 a = Cj+Cjj »+c 3 >' — w ' *"'" «" 3 + w * * co * «"' or 

a _ c,+c 4 »-|-c 3 v I wVK^Cy' sin ct + v 2 cos w) . _. 

* r:p (i +c v) (1-ho"*" ( ^ ^ 

Where Ci, c 2 , c 3 are constant quantities, of which two are arbitrary 
with respect to the given equation. 

The former fraction = c ' ~~P t + c ?-' — c 3 and therefore maybe written 
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1 + »* C 3 

The fraction (w) may be readily resolved into two fractions of 
the form 

P+? » + »+« 

where /> + r = o g + * c 2 = — m 2 Kc sin *■ 
? + s = o p + r c 8 =r — m 2 JC cos » 

1i — — - m> ^coso _^ m* Ac sin « 

^ ~~ 1— c* ' * ~" 1— c* 

„ nt'Keotv „ m* isfc sin «r 

r= , s = 



the fraction ^j" may be considered as contained in 

iqp-r * c x and c 2 being arbitrary. 

whence z = c -^r + TF?7 ~~ c * 
Therefore, when w is odd and » = o 

71—1 "— 1 

d»M=d» s= (—1) — c 2 d»" + (— 1) 2 9 C 1 - 1 d»" 

when n is even and v = o 

d» M = d n s =(—1)^ d » n + (— 1)> c* d i* 
These equations are the nth particular fluxions of 
w - c, cos » + c 2 sin » + p cos c » + ■£- sin c » — c 3 
comparing the second fluxion of this equation with the given 

equation 

— c} = G 
Therefore substituting for p and q their values above given 

u-c 1 cosv + c 2 sin »— £^. cos (c »— ») + G 

Case 2. When c = 1 

Equation (6) of this article gives 

h 2 
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a = TT'" — c - 5 — ( i+>') s i; 

let the term « be represented by 

let also the part of u corresponding to v be called w. 
1. n being even and," *= o 

2. n being odd and v = o 

n— 3 

d» v = (— \)~ (*T-)qdv n . Therefore d" w =(— I)^*") rf r"- 1 </» 



n— 1 

n— 1 



Or <£"U> = (—1) 2 ( 2 d l* n Cffrr-|(fl* n ) 

according as n is even or odd. 

These two equations are the particular divided fluxions of 

w= — - v sin v + i- v cos v — !L sin v 

2 ~ 2 2 

the latter term may be neglected, being contained in 

c 2 sin v. Hence substituting for p and g 

u = Ci cos v + c 2 sin »» — i m 2 -K k sin (c — ») + G 

Case 3. When c = 1 and m 2 is a very small quantity, the following 
integration will be exact to the lirst power of mz. 

By the same substitution as in the preceding case 

now c < _i>>' r ,/J p>* 1 _ _£, y T 

«i(i+» ) 
nearly = r+Tr+^yT 

Hence * = i + ^i +£ \ , + l+'/i + yp" nearly 



45 

Then d"u =<T z = c t (-iy O + p\* dv, n even and v = o 

n — 2 n — 1 

d*u = tfz = c, (_l/"§Vl +1 )~d / , » odd and 9 = o 
These two equations are the particular fluxions of 
m = c, cos (1 +£\» + c, sin (1 +i\ v + G 

where p ■=. m 1 Kcos «r and q — m' K sin *■ 

A first approximation therefore gives the value of « in period- 
ical terms, instead of in terms without the periodical signs as occur 
in the complete solution of the differential equation obtained in 
case the 2d. 

(IV.) 
Further application to the latnar Orbit. 

It was found by Equation (18) Art. (II.) that without the dis- 
turbing force of the Sun 
« = If 1 + e cos ( » — r)\ = I + _l cos » cos «■ + _£_ sin v sin *, re- 

<t V. J a a a 

garding only the first power of the excentricity. 

Also without the disturbing force of the Sun by the preceding 
article the integral of equation (d) of article (II.) becomes 
u = c , cos f + c t sin » + i_ 

a 

comparing these values of u 

c, = ± cos r c t = .£ sin w 

Hence in case 3 of the preceding article, because 

KS e 1,1 /. Sin' e cos a- 3 m z e sin «■ 
= — £- and therefore p = — s- » o = — — n 

2a r 2a * 2a 
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c, cos CI + I y ' + C sin " ~ )* ' 

= - costt cos(1- — )*»+-sinirsin(l— _J% = 

Consequently by the application of the first and third cases of 
the preceding article, the integration of equation (d) 

gives u = 7TTIl - 27 +- cos(/ 1 — -) »—■ ) + — cos (2 p-8iw) 
+ £l cos 0— 2 m » + *>+-7- « cos ( 3 *— 2 m » — *) (/} 

where ^ <0) ==-^(^) (T ^ z ^ und therefore is of the order 

of i«' 

^ = 3w'j 2+ W i-(i-2w)» and tnerefore on account of the 

divisor 1 — (1 — 2 m)* = 4 m— 4 m* is of the order of m 

^t( 2 > = 3m'| 2_2 m r l— la-law)* ant * therefore is of the order 

of m 1 

Hence u — — "-— cos (v — + — cos [r — m v + -*) will be a new 

value of w, exact to the order of m. 

It is clear by the third term of the value of u given by equa- 
tion (f) that when (1 — -~) * — *" = o or a multiple of the circum- 
ference, the moon is at perigee, not regarding the periodic terms 
depending on the place of the sun, and therefore the mean motion 

of the moon : mean motion of its perigee : : 1 : ~ The pro- 
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gressive motion of the perigee thus found is, as is well known, 
only about half its real motion. 

We must next investigate the effect of the new substitution 
u = — (- — cos (V— tt) + — cos (V — m * + *■) on the terms of the 

equation (c) of art. (1.) and thence deduce a new equation (rf) the 
integration of which will give a new value of «. For this 
purpose it will only be necessary to compute the variation of 

equation (d) arising from a variation i u = ^- cos (v— m * + w) 

and as the object of our enquiry is the mean motion of the Lunar 
perigee, it will be only necessary to compute the new terms of the 
equation (d) of the form B cos (v — »■), for from a similar term, viz. 

— ^^cos (»— *-) arose the term cos((l — — )*— *). The inte- 
gration given in case 3 of the preceding article will then at once 
give the coefficient from which the mean motion of the perigee 
is deduced, as far as the new approximation to the value of u is 
concerned. 

To compute the variation of the equation (d) as to the above- 
mentioned terms. 

» = m v — 2 e m sin {» — *) 
This gives 3 '» = — 2 e m sin (•—*•), supposing d u = — cos (»—*•) 

Hence nearly £ v = — 2 -4 e m sin (» — m v + v) 
when hu - — cos (» — 2 m v + *■) 
Substituting these values of £ u and £ >, 
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The first term or — i"lVl contains no term of the form 

2 a 3 h u* 

B cos (»— ?r) 

Drm rtl" ■ __ 



The second term or — ^rr* cos (2 » — 2m c) gives a term 



Sw'ejl 9m* A e 



4>a+h u* ~~~ *a 

The third term gives a term 



cos(» — cr) (I) 



ha*u 3 



c — it) = cos (t — x) (2) 

e^t. rfu g . n (2> _ 2 m f) __ w rfu gin 2 ,_ 2 

— a 4 u* tit v ' '2a 4 u* rf» v ' 

+ .!^gLa; coB(2»-2mO 

1 a*« 4 flv v 

The first and third terms do not afford any term of the required 
form B cos (■< — ») 

hxA — 5^' T, sra (2»— 2m>) gives 

i (1— 2m) ~ £ cos (r-r) (3) 

3.. The variation of 

$? + « }/& = -» |£? +.}/£*.*.(*-*•) 

= (substituting - instead of ■— -f u )jT ™ u " ^" sm (2» — 2») 

The first two terms give a term 

— (1 + m) A e cos (»— jt) (4) 

The third term gives a term 

- — • Srf-ir cos o-" ( 5 ) 
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Hence collecting the terms (1) (2) (3) (4) (5) the co-eflicient of 
the term e cos (« — w) in the variation of the equation (d) obtained 

by iu — — cos (^ — 2m»+j) will be found ■= 

57 (5+ 5 m) A e cos (k — *) 

consequently the coefficient of e cos (v — r) in the differential equa- 
tion arising from the substitution of 

u — -A — cos (v — t) + — cos (» — m v + x ) 

a ' a \ J a v J 

is thus found = — ^ (1 + (5+5m) A) 

Hence by case 3 Art. (Ill) the corresponding term in the in- 
tegral or value of u = e cos f (1 — (1+ (5+5 m) A)*—*) 

Therefore the mean motion of the perigee is 

i m 2 (1+ (5+ 5m) A) that of the Moon being unity. 

<< = * (c=) (T + ?£) = i«(5+19») nearly; 

Periodic time moon , ««*»»«<-> 

now m = p • ,v 7 it nearly = 0,0748. 

Periodic time earth J 

Therefore A — , 1801 nearly. 

Hence the mean motion of the perigee thus found by the se- 
cond approximation =, 00826 or 2°.58',4 each revolution, the ob- 
served mean motion = , 00845 = 3. 2, 5 each revolution : the dif- 
ference is about ~~\h part of the whole, whereas the first approxi- 
mation was only about one half of the whole. 

It is evident, that if we had also found the variations arising 
from the other terms of u besides Ae cos (r - m v+w), the results 
would, (as will easily appear by considering the values of 
these quantities in comparison of Ae cos {a — mv+ff) and the na- 
ture of the combinations required to form the angle (»-- v) have 
had only a small effect on the mean motion of the perigee. 

The above, it is conceived, is sufficient for two purposes. 

VOL. XIII. I 
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1. To shew by a compendious process that the mean motion of 
the Lunar perigee may be satisfactorily accounted for by the New- 
tonian Theory of Gravitation. 

2. To shew that the term e cos (»— w) occurring in the first ap- 
proximation for the value of u becomes by succeeding approxima- 
tions of the form e cos (c»— *) and therefore if we use this form in 
the first approximation and then add other terms of the proper 
combination of angles with indeterminate coefficients, as Laplace 
has done, considerable advantage will be thus afforded of obtaining a 
sufficiently exact value of u, for the purpose of an exact determination 
of the Lunar orbit. No farther difficulty afterwards occurs ex- 
cept that arising from the great length of the computations. 

It is remarkable that in the above expressions for the mean 
motion of the perigee, the part f m 2 arising from the first 
approximation is nearly equal to i m 2 (5 ±5m ) A ; because 
m = ,0748 and A =, 1814 and therefore (5 + 5m) A = 0,97 or 
unity nearly. But this coincidence is merely accidental, had the 
periodic time of the moon been double what it is, the first approxi- 
mation would not have given ;- of the whole motion, but had the 
Moon revolved about the Earth in one day, the Newtonian re- 
sult, according to the principles of the 9th sect, of the first book 
of the Principia, would have been within T *_ part of the whole, 
instead of being within one half only. 

Hence we see the results of the investigations of Machin, 
Walmsly, Mathew Stewart and Frisi are all quite erroneous, as 

their methods give the mean motion of the perigee = -^ nearly, 

which only affords an exact result when the periodic time of the 
Moon = T l T that of the Earth nearly. 
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The common error into which they have fallen is, as appears to 
me, pointed out in a paper of mine in the 8th volume of the trans- 
actions of this academy. 

Although the above investigation is only exact to the first power 
of the excentricity of the lunar orbit, yet it appears that the ex- 
centricity does not enter the above expression for the motion of 
the perigee, and therefore it follows that the excentricity of the 
lunar orbit has but little effect in the motion of the perigee. 
In fact only the second powers of the excentricities of the solar 
and lunar orbits and the second power of the tangent of the in- 
clination of the lunar orbit to the plane of the ecliptic enter 
into the more exact expression of the mean motion of the lunar 
perigee. 
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